Phenotypic variability in traits related to flight dispersal in the wing dimorphic species Triatoma guasayana

Background Triatoma guasayana is considered an emerging vector of Chagas disease in the Southern Cone of South America. The presence of a triatomine population with brachypterous individuals, in which both wings are reduced, has recently been reported for this species. The aim of the present study was to determine if flight-related traits varied across populations, if these traits could explain differences in flight capacity across populations and if flight-related traits are associated with geographic and/or climatic variation. Methods The study involved 66 male T. guasayana specimens from 10 triatomine populations. Digital images of wing, head and pronotum were used to estimate linear and geometric morphometric variables. Variations in size and shape were analysed using one-way analysis of variance and canonical variate analysis (CVA), respectively. Mantel tests were applied to analyse the relationship between morphometric and geographic distances, and the association between size measurements was analysed using Pearson’s correlation. We explored covariation between size and shape variables using partial least square analyses (PLS). The association of geographic and climatic variables with size measurements was tested using linear regression analyses. We performed PLS analyses for shape measurements. Results Wing size differed significantly across triatomine populations. The CVA showed that wing shape of the brachypterous population is well discriminated from that of the other populations. The Mantel test showed a positive and significant association between wing shape and geographic distances. The heads of the brachypterous population were significantly larger than those of the other populations. Similar to wing shape, the head shape of the brachypterous population was well discriminated from those of the other populations. Pronotum width did not show significant differences across populations. Geographic and climatic factors were associated with size and shape of both the wing and head, but not with pronotum width. Conclusions Most of the traits related to flight dispersal varied across populations. Wing shape and head shape were found to be better markers for differentiated morphological variation across populations. Head measurements also varied in accordance with this condition. Geographic and climatic variables were associated with most of the flight-related traits. Graphical Abstract Supplementary Information The online version contains supplementary material available at 10.1186/s13071-022-05570-7.

also varied in accordance with this condition. Geographic and climatic variables were associated with most of the flight-related traits.
Keywords: Chagas disease, Head, Pronotum, Wing, Linear morphometrics, Geometric morphometrics, Wing dimorphism Background In many organisms, environmental heterogeneity and ecological gradients can generate phenotypic variation [1]. Phenotypic patterns of morphological variations within insect species might reflect different dispersal patterns and habitat availability. Intraspecific morphometric analyses should focus on samples from a broad geographical range to understand variations in morphometry of individuals that occur in diverse environments. Additionally, to understand how environmental heterogeneity contributes to changes in phenotypic variation, knowledge on whether some morpho-functional mechanisms generate this pattern of variation is desirable [2].
The subfamily Triatominae (Hemiptera: Reduviidae) includes hematophagous insects that act as vectors of Chagas disease, one of the most important parasitic diseases in Latin America [3]. Currently, there are 157 living species [4] of triatomines, which differ in their epidemiological importance and in aspects of their biology [5]. In recent years, several studies have reported the occurrence of sylvatic triatomines dispersing actively to domestic environments in Argentina, Brazil and Paraguay [6][7][8].
Although sylvatic species are considered to be of secondary epidemiological importance, they are responsible for maintaining the sylvatic transmission cycle of Trypanosoma cruzi [9].
Triatoma guasayana is a sylvatic species that invades peridomestic and domestic habitats. It is distributed in the dry western Chaco region of Argentina, Bolivia and Paraguay. Some laboratory and field studies have shown that this triatomine species has vectorial capacity [8,[10][11][12][13]. In sylvatic ecotopes, T. guasayana has been found in dry cacti, bromeliads and fallen dry logs, and in peridomestic habitats, it is associated with chicken coops, goat or sheep corrals, piled materials and orchard fences [14][15][16]. Triatoma guasayana is traditionally considered to be part of the "sordida subcomplex" [17]. However, phylogenetic analysis indicated that T. guasayana is more closely related to the "rubrovaria subcomplex" than to members of the "sordida subcomplex" [18,19]. According to Gorla and Noireau [20], T. guasayana, among other species, may be considered an emerging vector in the Southern Cone of South America.
It has been postulated that flight dispersal is the main mechanism of invasion, infestation and/or recolonization of houses by adult triatomines [21,22]. Nutritional and reproductive status as well as population density and environmental conditions are known factors that modulate flight dispersal in member of the Triatominae [23][24][25][26]. There are interspecific differences in terms of dispersal capacity by flight; for example, the number of flying individuals was found to be higher in T. guasayana than in T. infestans and other secondary triatomines [12]. Flight dispersal may vary markedly even within the same species. Some species present polymorphisms in wing length and/or in the development of the flight muscles that affect the active dispersal capacity [27,28]. Wing polymorphism implies variation in the size of the wings or in their presence/absence [27]. The diversity of wing morphotypes is usually associated with other morphological attributes, such as head shape, thorax size and olfactory capacity. In Mepraia spinolai, a species with wing dimorphism, macropterous individuals (with normal wing size) exhibited a greater development of the thorax than micropterous specimens [29]. For T. guasayana, the presence of a population with brachypterous individuals, in which both wings are reduced, has recently been reported [30]. In this study, changes in the shape of the head according to wing characteristics were also determined, with macropterous individuals found to have a shorter head and a greater distance between eyes [30]. These observations in M. spinolai and T. guasayana suggest that the varieties of wing morphotypes imply a diversity of body structures. Polymorphism in flight muscles has been reported and includes developmental changes that affect dispersal capacity [31]. In some triatomine species, a relationship between wing size and thorax development has been demonstrated. Hybrids between Triatoma sherlocki and Triatoma juazeirensis have an increased dispersal capacity that coincides with a wider thorax, despite having wings of intermediate size.
A wider thorax may imply greater development of the wing muscles. Mepraia spinolai showed an increase in thorax size as a function of wing development [29].
Taking into account the existence of wing dimorphism in T. guasayana and that this condition determines flight functionality, the aim of this study was to understand: (i) if traits related to dispersal characteristics varied across populations; (ii) if size and/or shape variation of flightrelated traits could explain differences in flight capacity across populations; and (iii) if flight-related traits are associated with geographic and/or climatic variation. Gigena et al. Parasites & Vectors (2023) 16:8

Insects
The study involved 66 males from 10 populations of T. guasayana (Table 1; Fig. 1). Only males were studied because the sample included more males than females, and because other studies suggest sexual dimorphism ( [30] and J Nattero, unpublished data). For the Cochabamba population, from the 16 collected males, wings of five individuals could be included since the remaining individuals had not been properly preserved. All populations included in this study were collected from peridomiciliary structures. All adults from the Cochabamba population are brachypterous, a condition first reported by Hernández et al. [30] (Fig. 2a). For this population, the dispersal capacity in the field was evaluated with light traps, active searching and live-bait mouse traps (also known as Noireau traps). Although individuals were captured by active searching and with Norieau traps, no

Geographic and climatic characterization
To obtain the geographic location for each population, we recorded latitude, longitude and altitude. The geographical coordinates of the 10 populations were imported into the Google Earth Engine platform to extract the climatic variables for each site. We used ERA5 (Copernicus Climate Data Store; https:// cds. clima te. coper nicus. eu/# !/ home) to obtain the climate data with a spatial resolution of 11 × 11 km. Of the 50 climate variables available, in this work we only included data from annual rainfall and annual mean temperature. We also considered the average monthly dew point temperature and the average monthly temperature and estimated the relative humidity (RH), including these data as a third variable to characterize climate variation. RH was calculated as [32]: where RH indicates relative humidity, T represents the average monthly temperature and Td is the average monthly dew point temperature. We only included these three environmental variables because they are determinants of flight dispersal in insects [33,34].

Linear morphometry
Linear measurements of the head and pronotum were made using the free software package tpsDig2 version 2.31 (http:// life. bio. sunysb. edu/ morph/). The linear measurements of the heads included anteocular distance (AD), width between antenniferous tubercles (WAT), preocular width (PW) and maximum width across the eyes (EW) (Fig. 2b). For the pronotum, we measured the maximum thorax width (MTW) (Fig. 2c). All linear measurements showed a normal distribution.

Geometric morphometry
A landmark-based approach was applied to study geometric morphometrics of the heads and right forewings. We defined and collected five coplanar type II landmarks of the ventral view of the head, and eight type I landmarks of the wings using the software package tpsdig2 (Fig. 2d, e). We extracted data on head and wing shape with a generalized full Procrustes fit and a projection to shape tangent space [35]. We used Procrustes coordinates as shape variables. Centroid size (CS; i.e. the square root of the sum of squared distances from each landmark to the centroid of the configuration) was computed as a measure of wing size and head. All of these steps of the morphometric analysis were performed using the free software MorphoJ version 1.07a [36].

Statistical analysis
Linear and geometric variables of the head, pronotum and wings were analysed separately. Variations across populations for the single linear measurements of the wing, head and pronotum were explored via one-way analysis of variance (ANOVA) and Tukey's post-hoc tests for each variable, using the software InfoStat version 2016 [37]. For shape measurements of the head and wings, we calculated Procrustes distances between pairs of populations and evaluated their significance via a non-parametric test based on permutations (1000 runs), using MorphoJ [36]. We then performed canonical variate analyses (CVA) and calculated the percentage of phenotypic similarity between pairs of populations using the cross-check test of discriminant analysis [38] in InfoStat software [37]. The Procrustes distances were represented in unrooted neighbour-joining (NJ) trees using the free software MEGA X 10.2.6 [39]. We calculated bootstrap values for 1000 replications following Ascarranuz et al. [40]. The relationship of geographic distances with linear and geometric morphometric measurements of head, pronotum and wing was analysed using Mantel tests with the free software PASSAGE 2 2.0.11.6 [41]. To construct the morphological distance matrices between populations, we used Procrustes distances for shape measurements and Euclidean distances for size measurements.
The association between wing, head and pronotum size was explored via a Pearson's correlation, using InfoStat software [37]. We analysed the covariation between wing, head and pronotum shape and between size and shape using partial least square analyses (PLS) in MorphoJ [36].
To summarize the variation and describe the environmental characterization of the six geographic and climatic factors (i.e. latitude, longitude, altitude, annual rainfall, annual mean temperature and RH), we used a principal component analysis (PCA) with R package FactorMineR 2.4 [43]. The first two principal component (PC) scores described 95.04% of the total variance (59.72-35.69% for PC1 and PC2, respectively). In PC1, the factors with relevant coefficients (> 0.7) were latitude, altitude, temperature, precipitation and RH (Additional file 1: Table S1). In PC2, only longitude was an important factor (Additional file 1: Table S1). We explored the association between PC1 and PC2 and measurements of the head, wing and pronotum size across the 10 populations studied using a multiple linear regression analysis [43]. The covariation between geographic and climate variables and wing and head shape variation was explored using PLS.

Changes in wing size across populations
Significant differences in wing size were observed among populations (F (9, 42) = 2.99, P = 0.008) (Fig. 3). Tukey's post-hoc tests showed significant differences between the Loreto population and two other populations: Cochabamba, the brachypterous population, and Independencia (P < 0.05) (Fig. 3). The Loreto population presented a smaller wing size than the Cochabamba and Independencia populations, whereas the latter two had the largest wings of the 10 populations. The Mantel test showed a negative and significant association between wing size and geographic distances (Z = 73,661.06, r = − 0.05, P = 0.04), indicating that the larger the differences in wing size, the smaller the geographic distance.

Wing shape of the brachypterous population differs from that of the other populations
Wing shape differed among populations. The first two discriminant factors of CVA (CV1, CV2) explained 68.13% of the total variation (CV1: 42.60, CV2: 25.53%). The brachypterous population, Cochabamba, was found to be well discriminated from the other populations in the space of the first discriminant factor (CV1) (Fig. 4a). All individuals from this population, as well as those of the Campo Alegre, Cruz del Eje, Ischilín and Valle Fértil populations, were correctly assigned (0% of misclassified individuals) (Additional file 1: Table S2). In contrast, the San Martín population exhibited 60% of misclassified individuals (Additional file 1: Table S2). Procrustes distances between populations were significantly different (P < 0.05) for all of the comparisons that included the Cruz del Eje and Cochabamba populations (Additional file 1: Table S3). The NJ tree based on Procrustes distances is presented in Fig. 5a, which shows that the only bootstrap-supported clade was the one that included the Independencia and Cruz del Eje populations and the Ischilín and Campo Alegre populations (Fig. 5a). The Mantel test showed a positive and significant association between wing shape and geographic distances (Z = 2483.882, r = 0.44, P = 0.001), indicating that the smaller the differences in wing shape, the smaller the geographic distances.

Head size of the brachypterous population is larger than that of the other populations
The ANOVAs performed for the four head linear measurements showed significant differences in all cases ( Table 2). Tukey's post-hoc tests showed that the brachypterous population had the highest values for all four measurements, although these were significantly different from those of the other nine populations only for the anteocular distance (P < 0.05) ( Table 2). In the Mantel test, head size showed a positive and significant correlation with geographic distances (Z = 83,069.24, r = 0.324, P = 0.001,) indicating that the smaller the differences in head size, the smaller the geographic distance.

Head shape of the brachypterous population is different from that of the other populations
The first two discriminant factors accumulated 96.72% of the total variation (90.18 and 6.62% for the first and second axes, respectively). The space delimited by the first   Table S2), while The population from 12 de Junio exhibited the highest number of misclassified individuals (85.71%) (Additional file 1: Table S2). Procrustes' distances between populations were significantly different (P < 0.05) for all pairs of populations that included Cochabamba (Additional file 1: Table S3). The NJ tree based on Procrustes distances is presented in Fig. 5b. Bootstrap values supported different groups from different geographic areas; in particular, the brachypterous population was in a group that was separated from the rest (Fig. 5b). The Mantel test showed a positive and significant association between head shape and geographic distances (Z = 4499.254, r = 0.19, P = 0.005), indicating that the smaller the differences in head shape, the smaller the geographic distance.

Pronotum width is similar across populations
The ANOVA for the pronotum width did not show significant differences across populations (F (9, 51) = 1.78, P = 0.094). The Mantel test showed a negative and significant correlation with geographic distances (Z = 71,936.392, r = − 0.147, P = 0.024), indicating that the larger the differences in pronotum width, the smaller the geographic distance.

Wing, head and pronotum measurements are associated
Pearson's correlations showed positive and significant associations between wing, head and pronotum size measurements ( Table 3). The PLS analyses showed no covariation between wing and head shape, wing shape and head size, or between wing shape and pronotum width (Table 3). Head shape showed significant covariation with pronotum width and wing size (Table 3).

Linear and shape measurements are associated with geographic and climatic variables
Multiple linear regression analysis showed that head linear measurements were significantly associated with geographic and environmental factors with both PC1 and PC2 coefficient (with the exception of width between  Table 3 Results of Pearson's correlation between linear measurements of wing, head and pronotum, and partial least square analyses between shape measurements of wing and head and between size and shape measurements of wings, head and pronotum in T.  (Table 4). In contrast, pronotum width and wing size were not associated with geographic or climatic factors (Table 4). PLS analysis showed significant covariation of head and wing shape with all geographic and climatic variables (Table 5).

Discussion
Triatoma guasayana exhibits a good dispersal capacity. In a study comparing this species with other triatominae species (T. infestans, T. platensis and T. eratyrusiformis) from the arid Chaco of Argentina, T. guasayana was found to have the highest average number of flying individuals [12]. Adult T. guasayana bugs may invade a wide variety of ecotopes and have a high potential for night flight dispersal during the hot dry months (August-November) [21,44]. Since dispersal is linked to flight, the existence of wing dimorphism at the population level provides a suitable model to study not only evolutionary patterns, such as dispersal processes and population structure [45,46], but also a range of ecological factors, including habitat structure and climate [47,48].
Our results showed differences in wing size across T. guasayana populations. Although all triatomine males collected from Cochabamba showed a brachypterous condition, individuals from this population as well as those of the Independencia population exhibited a bigger wing size than the macropterous populations.
Individuals from the brachypterous population had significantly greater body size than individuals from the two macropterous populations of Ischilín and Cruz del Eje (F (2,11) = 6.76, P = 0.012). This result suggests that the individuals from the Cochabamba population are larger and that their body is comparatively larger than their wings, showing the presence of brachyptery in this population. The same tendency has been observed in other insects, such as in carabids, with brachypterous individuals exhibiting a larger body size than their macropterous counterparts. This variation is attributed to differences in dispersal strategies, taking into account that brachypterous individuals disperse by running [49].
Regarding wing shape, our results showed differentiation between the brachypterous population and the other nine populations included in this study. This result suggests that differences in wing shape are associated with the brachyptery condition, and that this pattern could be influenced by the differences in flight capacity across populations. Differences in the shape but not in the size of a structure are considered to be determined by genetic processes [50]. Our results also showed that there is no covariation between wing shape and size measurements of the wing, head and pronotum. These results suggest a lack of functional integration (i.e. coordination among parts) of wing shape and size measurements across populations, and may be related to wing dimorphism. Studies combining biomechanics of flight and morphometric   16:8 approaches could contribute to a better understanding of the functional implications of morphological changes associated with wing dimorphism [61].
Insects use a combination of sensory inputs, including cephalic sensilla, ocelli, antennal sensilla and eyes, to control flight activity [52][53][54]. Our results show that, in general terms, the brachypterous population presents a larger head size than the remaining populations. For T. guasayana, anterocular distance and distance between eyes have been suggested as the differences in the head that could be associated with flight capacity; it has also been proposed that shorter heads are related to changes in the position of the ocelli [30]. We found significantly longer head measurements for the brachypterous population, which agrees with the association suggested by [30]. Regarding head shape, our results showed that the head shape of the brachypterous population differed from the that of the other nine populations. The NJ tree for Procrustes distances showed similarities between populations that are geographically close.
Flight muscles are located in the pronotum. Our results showed that pronotum width was the only measured trait that did not differ between populations. There are alternative, non-mutually exclusive explanations that might account for this result. First, among insects, there is muscle polymorphism associated with the degree of muscle development. A greater pronotum width does not necessarily indicate a greater development of these muscles, as has been reported for other insect species [55]. Alary polymorphism is not necessarily accompanied with flight muscle development [30]. Second, there are functional or developmental constraints that might limit the potential for evolutionary change of pronotum width and reduce the potential of the population to respond to selection or to undergo non-adaptive evolution by drift [56].
The association of morphological distance with geographic distances between populations showed contrasting results. While wing shape and head size and shape exhibited positive associations, wing and pronotum size showed negative associations. Thus, these results might imply geographic separation as a key factor in morphological differentiation of genetic exchange between populations. A positive association would imply that geographically closer populations are more similar to each other. It would also indicate morphological spatial structuring for wings and heads between the populations of this species. This positive spatial association would adjust to a model of isolation by distance, in which differentiation increases with increasing distance. By contrast, the sizes of the wing and pronotum showed negative associations, implying that the greater the distance between populations, the smaller the difference in wing and pronotum size, respectively. These results should be interpreted with caution. In some cases, the correlation values were low, such as those for wing size and shape and for pronotum width.
The association of geographic and climatic variables with size and shape measurements showed that the former variables better explain head size, head shape and wing shape than the other measurements. For wing size and pronotum width, only longitude showed association. This means that, considering geographic variables, all morphological measures followed an east-west variation, and that head size, head shape and wing shape additionally followed a north-south variation and an altitudinal variation. Considering the broad-scale geographic variation of the sampled populations, variation in geographic variables also determines variation in climatic conditions. Wing and head shape exhibit significant phenotypic plasticity in response to developmental temperature, and this plastic response (hot-to-cold) mirrors a pattern of geographic differentiation, as might be expected under a synergistic scenario of adaptive plasticity and adaptive genetic divergence [57,58]. Similar results were found in the yellow dung fly Scathophaga stercoraria (Diptera: Scathophagidae) from the New World, with clines of variation associated with wing shape but not with wing size [59]. There are also examples in Drosophila flies of no association between wing size and temperature variation (13 °C and 22 °C) [60].
In M. spinolai (sub Triatoma spinolai [61]) collected from various sites in northern Chile, adults from coastal populations were always wingless, while inland populations showed both wing and wingless adults. Winged males would be able to disperse, although this event does not guarantee that fliers will land in favourable sites for feeding and/or mating [61]. Brachyptery in T. guasayana is a heritable trait for at least three laboratory generations [35]. In other Hemiptera, the genetic base of brachytherapy has been found to be polygenic [62]. The morphological characteristics of the heads and wings observed in the Cochabamba population suggest a lack of dispersal by flight; therefore, the population would be genetically isolated. The differentiation observed in all of the characters would support this assumption, especially for those characters related to shape variation, which suggest genetic differentiation and, as mentioned above, this genetic differentiation between populations could be indicating different modes of dispersion. On the other hand, our results show that both head shape and wing shape correlate with geographic and climatic variables. Other insect species, such as Drosophila subobscura, or the fly Scathophaga stercoraria, also show a geographic pattern of morphological differentiation, including significant clinal variation in wing shape [52,60].